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Abstract 

Buckminsterfullerene compounds exibit remarkable physics at low 
temperatures, e.g. high temperature superconductivity in alkali-full- 
-erenes, and ferromagnetism in TDAE- Cgo- Here we review recent 
theoretical studies of electron correlations in these compounds. In 
particular, we discuss models of electron-vibron interactions, electron- 
electron interactions, and intermolecular hopping. We show that the 
origin of novel electronic phases lies in local degeneracies of Ceo; ^ 
direct consequence of the high molecular symmetry. 



1 Introduction 

The synthesis of buckminsterfullerene ( Cgo) into molecular crystals with elec- 
tron donors has resulted in materials with surprising electronic properties]^]. 

Ceo is a truncated icosahedron. From a physicist's standpoint, the charged 
molecule is fundamentally interesting, because the high molecular symmetry 
gives rise to degeneracies in both electronic and vibrational systems. Thus, 
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the molecule is very sensitive to perturbations. In particular, electron- 
phonon and electron-electron interactions are expected to produce highly 
correlated ground states and excitations. 

First, superconductivity was discovered in alkali-fullerenes, A3C60, (A=K, 
Cs, Rb) at relatively high temperatures (Tc< 33°K P]), compared, say, to 
intercalated graphite. Soon thereafter ferromagnetism was found in TDAE- 
Ceo at Tc~ 16°K P, ||, where TDAE is tetrakis(dimethylamino)ethylene, 
C2N4( 0133)3 (the stoichiometry is a 1:1 ratio of TDAE to Ceo)- 

Some striking aspects of TDAE- Ceo are its relatively large value of Tc- for 
an organic ferromagnet - and and its nonmetallic conductivity, suggestive of 
Mott-Hubbard localization[0. Ferromagnetism is thus particularly unusual, 
since superexchange naively predicts antiferromagnetic interactions between 
localized spins in many such systems. 

An obvious question is raised: Are superconductivity in A^Cqq and ferro- 
magnetism in TDAE- Ceo related? While superconductivity involves effec- 
tively attractive interactions, magnetism is usually believed to result from 
repulsive Coulomb forces.^ The similarity between the materials is that both 
involve partially filled conduction bands made of ti^-orbitals of Ceo- Also 
one expects similar intramolecular electron-vibron and electron-electron in- 
teractions. 

The primary differences between the two systems are in their crystalline 
symmetries and Ceo ionizations. While A3C60 is an FCC crystal with cubic 
symmetry in which the Ceo molecules are triply ionized (Ceo ), TDAE- Ceo 
has a c-centered monoclinic unit cell0, which gives rise to preferred hopping 
along the c-axis, and singly ionized Ceo- 

Here we review recent work on the local molecular interactions which are 

separated into two major contributions: the electron-electron pseudopoten- 

tials (Section ^ , and the electron-vibron interactions (Section The latter 

""^A similar mystery underlies the proximity of high TcSuperconductivity to antiferro- 
magnetism in the high- TcCuprates. 
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exhibit, in the semiclassical hmit, interesting dynamical Jahn Teller effects 
accompanied by non trivial Berry phases p, ^ |10] . We show that this effect 
enhances the attractive interaction due to electron-vibron interactions by a 
significant factor. This analysis, though admittedly qualitative, points at a 
plausible cause for the relatively high Tc's of A^Cqq. 

In Section ^ we discuss TDAE- Ceo- We derive an effective model for 
the low excitations of TDAE- Ceo using a muticomponent superexchange 



expansion about the Mott-Hubbard insulating phase [11 1. As a result, we 



obtain a rich phase diagram which includes spin ferromagnetism, spin density 
waves, and orbital ferromagnetism depending on the ratios of interaction 
pseudopotentials. 



2 The electron-electron pseudopotentials 

We consider a system of vr-electrons hopping on the truncated icosahedron 
(the soccer-ball) lattice of Ceo- The single particle Hamiltonian is 

^° =E^aE E (1) 

A ^J. s=u 

A/x denote the irreducible representations of the icosahedral group x\fj,{^) 
where f2 is a unit vector pointing at one of the soccer ball vertices. At low 

XXfii^) are closely related to the angular momentum functions Yi^i^)- 
However, the structure of the states near the Fermi level of 60 electrons 
is different. The soccer-ball lattice splits the degenerate spherical / = 5 
multiplet, and the Fermi level lies in the semiconducting gap between the 
5-fold and 3-fold degenerate A = hg,tiu orbitals respectively. 

The Coulomb interactions on a metallic spherical shell is parametrized 
by two dimensionless variables g,a. 

H^-^ = \ j dVLidVL2V{VLi-(l2) : p{(li)p{(l2) : 
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V = g 27rr™„ ^ — — - YI,,{(1,)Ylm{(1^) , (2) 

P = Es^lV's is the density operator where i/jI = Ea^ Xa/.(^)4^s- 9 = 
e^/(ei?eo) is the strength of the screened interaction, and R is the radius of 
the baU. The power law tail of V{r) depends on a such that as a decreases, 
the potential has shorter range. < a < 1 interpolates smoothly between a 
S function and an unscreened 1/r Coulomb potential: 

l^-=° =g27^r^,Jin^-n2) ; = ^ (3) 

' '2 I 

When a changes from to 1, the averaged interaction over an area of radius 
rmin (analogous to Hubbard's U) is held fixed. 

The interaction in the basis of Hq can thus be explicitly given in the 
second quantized form as 

2 V 



L 



.(2L + i; 



V rL A A' fiLX A' t „t „_ 

where summation of repeated indices is assumed. 

For a molecule with N electrons in a closed shell (above which there is 
a gap larger than the interaction strength), the zero frequency interaction 
between two electrons in the open shell is given by 

uxfi = E'^j^2 " 2i?Ar+i + Em (5) 

which we call "pseudopotentials" . A/i denote the representation of the + 2 
electron state. The computation of ma^ for the Hamiltonian + W^'^ 
is a difficult task which was undertaken using second order perturbation 
theory p!^, ^ and later using perturbative renormalization group up to three 



loops order 
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For the real compounds however, estimation of the true pseudopotentials 
requires incorporating the electron-vibron interactions (see the following Sec- 
tion ^) as well as crystal fields, screening from neighboring molecules and 
polarization effects of a electrons. As we shall see, the electronic properties 
of the conduction bands in the solid are very sensitive to the combination of 
electron-electron and electron-vibron pseudopotentials. Thus even ab-initio 
calculations (which involve certain approximations for the higher order 
effects of interactions) may not be accurate enough for precise predictions of 
the ground state phase diagram. Here we shall therefore look for qualitative 
effects and leave the electron-electron pseudopotentials as free parameters 
within an acceptable range of magnitudes. 



3 The electron-vibron problem and supercon- 
ductivity 

In this section we solve the problem of a single Hg vibron coupled to ti^ 
electrons in a C^^ molecule. The model is too simplified for quantitative 
predictions for Ceo, but it contains interesting novel physics which will be 
important for further studies of this system. 

Semiclassically, a dynamical Jahn-Teller effect occurs[|l^, |T^. For n = 
1,2,4,5, the molecule distorts unimodally, giving rise to a pseudo-angular 
momentum spectrum, plus three harmonic oscillators. For n = 3, there is a 
bimodal distortion, which generates a spectrum of a symmetric top rotator, 
plus two harmonic oscillators. The pseudo-rotations are subject to non trivial 



Berry phase effects |ll9[] , which determine the pseudo-angular momenta L, and 
thus the degeneracies and level ordering of the low lying states. Strong Berry 
phase effects seem to survive even at moderate and weak coupling as shown 
by the exact diagonalization results. 

We find at weak coupling that the pair binding energy is a factor of 
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5/2 larger than the classical JT effect, and a factor of three larger than the 
pairing interaction of Migdal-Eliashberg theory of superconductivity. This 
enhancement can be interpreted semiclassically as due to large zero point 
energy reduction of the pseudo-rotations. From the weak coupling point 
of view, this effect is due to degeneracies in both electronic and vibronic 
systems. 

Migdal's approximation neglects vertex corrections in the resummation 
of two-particle ladder diagrams. This is justified only in the retarded limit 
uo « tp- Here we have considered the opposite limit, where the molecular 
ground state energies are solved first, assuming that the JT relaxation time 
is of the same order, or faster than the inter molecular hopping time. In this 
regime, we have found therefore that Migdal's approximation substantially 
underestimates the pairing interaction, and Tc, for these ideal molecular solids 
[ pO[] . This large effect suggests that some of the enhancement is likely to carry 
over to the real case of A3C60 metals, where electron hopping t and vibron 
frequencies are of similar strength. 

3.1 The Model 

We consider a single Hg (five dimensional) vibrational multiplet which cou- 
ples to n = 1, ... 5 electrons in an open ti^ shell, tiu and Hg are the 
icosahedral group counterparts of the spherical harmonics {Yim}m=-i) 
{^2m}m=-2 respectively. By replacing the truncated icosahedron (soccer 
ball) symmetry group by the spherical group, we ignore lattice corrugation 
effects. These are expected to be small since they do not lift the degeneracies 
of the L = 1, 2 representations. 

The Hamiltonian is thus defined as||lU| 

H =H° + H^-" , (6) 
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where, 

= nujY,(bijbM + l) +ie- fi)J2cl,Cms . (7) 



M rns 



creates a vibron with azimuthal quantum number M, and cj^^ creates an 
electron of spin s in an orbital Yim- By setting /i — e we can discard the 
second term. 

The Hg vibration field is 

u{d) = -^{Y^^ml + Y2Min)bM) , (8) 
where is a unit vector on the sphere. The ti„ electron field is 

= Y: YU^)cl, . (9) 

m=—l 

The electron- vibron interaction is local and rotationally invariant. Its form is 
completely determined (up to an overall coupling constant g) by symmetry: 

H^-" cx g fdQu{Q)Y.tPl{Q)^s{^) ■ (10) 

s 

Using the relation 

JdCl YLMmYi^,{Cl)YirnM oc {-1)^'{L,-M\lmi;lm2) , (11) 

where (■ ■ ■) is a Clebsch-Gordan coefficient |2T|], yields the second quantized 
Hamiltonian 

s,M,m 

X (2, M|l, -m; 1, M + m)cl,CM+ms • (12) 

The coupling constant g is fixed by the convention of O'Brien, who studied 
first this kind of dynamical JT problem |2^. Representation (|12D is conve- 



nient for setting up an exact diagonalization program in the truncated Fock 
space. 



7 



3.2 The Real Representation 

The semiclassical expansion is simpler to derive in the real coordinates rep- 
resentation. The vibron coordinates are 

6 2 



?M = ^ E M,m{bl+i-irb.m) , (13) 



where 



,m^o — (2 sign(/x)) 2 {Sf^^rn + sign(/i)5^^_^) , 
M,,o = 5^,0. (14) 
{g^i} are coefficients of the real spherical functions 

V O m 

^i?e(F2H(i^)) /x=l,2 
^F2o(^^) /x = . (15) 

^/m(F2H(l^)) /x = -l,-2 
We also choose a real representation for the electrons 



V2 



cL = 4s ■ (16) 

Thus the Hamiltonian in the real representation is given by 



' -\/3g-2 go - \/3g2 -\/3g-i 



2 . 



V3gi A/3g-i -2go 



/c 
(17) 



8 



This form of the JT hamiltonian is well known [T^ . Since the Hamilto- 
nian is rotationally invariant, its eigenvalues are invariant under simultaneous 
0(3) rotations of the electronic and vibronic representations. 



3.3 Jahn— Teller Distortions 

In the classical limit, one can ignore the vibron derivative terms in (|1^, and 
treat q = {q^} as frozen coordinates in H^~'" . The coupling matrix in H'^~'" 
is diagonalized by [p^ : 



w] 



(z-^r 

2 + VSr 





T{w) 



(18) 



where 
T = 



cos tl) sin tl) ' 
— sin i\) cos i\} 
1 



/ COS0 sin^^ 

1 
\ sin Q cos Q 



coscp 


sm 


— sin (j) 


cos 














(19) 



w = {(f), 6, ip) are the three Euler angles of the 0(3) rotation matrix T. In 
the diagonal basis of (^), the electron energies depend only on two vibron 
coordinates: 





z 


\o/ 



(20) 



By rotating the vibron coordinates q to the diagonal basis using the L = 2 
rotation matrix D^"^^ [^], one obtains 



q^{r,z,zu) 



E 



M,,^D^i]^,{w)M-}^,q,iO) 



(21) 



where M^_^ was defined in (j^. 
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By (pTD, and the unitarity of D and M, is invariant under rotations 
of w. Thus, the adiabatic potential energy V depends only on r, and the 
occupation numbers of the electronic eigenstates n^, where Y^i^i = n. 

V{z,r, [ui]) = ^{z^ + r^) + {ni{z - VSr) + n2{z + A/3r) - n^2z^ . 

(22) 

V is minimized at the JT distortions (^n, ^n, at which the classical energy 
is given by 

=min V^(^„,f„,n,). (23) 

The JT distortions at different fillings are given in Table I. We define 0, 9 
as the longitude and latitude with respect to the diagonal frame ( "principal 
axes") labelled (1,2,3) (3 is at the north pole). z,f parametrize the Jahn- 
Teller distortion in the real representation (||), as 

{u^^{(l)) =|(3cos2^-l) +^sin2^cos(20). (24) 

In Table I we present the values of the ground state JT distortions at all 
electron fillings. We see that electron fillings n = 1, 2, 4, 5 have unimodal 
distortions which are symmetric about the 3 axis, while n = 3 has a bimodal, 
about the 3 and 1 axes. The two types of distortions are portrayed in Fig. 
1^. We depict the distortions of ( plf ) for the unimodal and bimodal cases. 

3.4 Semiclassical Quantization 

At finite coupling constant g, quantum fluctuations about the frozen JT 
distortion must be included. In order to carry out the semiclassical quan- 
tization, we define a natural set of five dimensional coordinates r, z, w. w 
parametrize the motion in the JT manifold (the valley in the "mexican hat" 
potential V) and r, z are transverse to the JT manifold, since V depends on 
them explicitly. 
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For the unimodal cases the kinetic energy is given by 

il^'~i(i' + r' + r2(2V')' + {e^ + smH<i)^)) . (25) 
and the semiclassical Hamiltonian is thus 

3 

H"" = /ia;X^(ata^ + i) , (26) 

7=1 

where L is an angular momentum operator, and are the three harmonic 
oscillator modes of r. The energies are given by 

E--' = huj |^^L(L + 1) + j:{n, + i) j . (27) 

The rotational part of the eigenfunctions is 

ntiQ) = yLmm im)^, (28) 

where Cl — {9,(j)) is a unit vector, and [[jT-isDQ is the electronic adiabatic 
ground state. It is a Fock state in the principal axes basis. In terms of the 
stationary Fock basis Ifriay]) where a — x,y, z, the adiabatic ground state is 

[rias] 

Each overlap is a Slater determinant which is a sum of n products of spherical 
harmonics 

(K.]IK])o =j:q,]Y,,,{n)Y,,M ■ --y.^M , (30) 

where C\i,] are constants. 
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Now we discuss how boundary conditions determine the allowed values 
of L. A reflection on the JT manifold is given by 

(l -tl. (31) 

Spherical harmonics are known to transform under reflection as 

YLm^{-l)''YLn.. (32) 

Thus, by (^) and (0), the electronic part of the wave function transforms 

as 

\Mn^i-mK])^n. (33) 

The reflection ( pi]) can be performed by moving on a continuous path on 
the sphere from any point to its opposite. It is easy to verify, using (^Tj) that 
this path is a closed orbit of q & R^: 

m qi-n) = m . (34) 

Thus we find that the electronic wave function yields a Berry phase factor of 
(— 1)" for rotations between opposite points on the sphere which correspond 
to closed orbits of q. In order to satisfy (^8])) using the invariance of the 
left hand side under reflection, the pseudorotational Yim wavefunction must 
cancel the electronic Berry phase. This amounts to a selection rule on L: 

(-1)^+" = 1 . (35) 

Thus, the ground state for n = 1 and 5 electrons has pseudo-angular mo- 
mentum L = 1 and finite zero point energy due to the non trivial Berry 
phases. 

3.4.1 Bimodal Distortion 

The analysis of the bimodal distortions n = 3 proceeds along similar lines. 
The quantization of the pseudo-rotational part is the quantum symmetric 
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top Hamiltonian. Fortunately, it is a well-known textbook problem (see e.g. 
Ref. 0)- The eigenfunctions of a rigid body rotator are the rotational 



matrices 

Dl'lM, (36) 

where L,m,k are quantum numbers of the commuting operators L'^,L^,L^ 
respectively. and are defined with respect to the fixed z axis and the 
co-rotating 1 axis respectively. The quantum numbers are in the ranges 

L = 0, 1, . . . oo 
m,k = -L, -L + 1, . . . L . (37) 

The remaining coordinates are two massive harmonic oscillators modes 

r = {r — f,z — z). (38) 

The semiclassical Hamiltonian is thus 

Tjrot f2 1 2 



2 

= nujY.{a\a^ + \), (39) 

7=1 



and its eigenvalues are 

E'' = (^-^.HL + 1) - ^fc^ + |:(n, + i) j . (40) 

The rotational eigenfunctions are explicitly dependent on w as 

^ltM = Dlil{w)YlM^- (41) 

is 

Unlike the unimodal case, in the bimodal case no single reflection fully 
classifies the symmetry of the wavefunction. However, one can obtain definite 
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sign factors by transporting the electronic ground state in certain orbits. We 
define the rotations of tt about principle axis U as Cj. The Berry phases 
associated with these rotations can be read directly from the rotation matrix 
T in Eq. (|T9|) . For example: ioi ip^ip + (C3), the states |1) and |2) get 
multiplied by (—1). 

Since D^^\ transform as Y^k under Cj, it is easy to determine the sign 
factors of the pseudorotational wavef unction. The results are given below: 

C, : |1,0,2)^ |1,0,2)^, Ci : ^i^^ ^ (-l)'l^S 



C2 : |1,0,2)„ ^ -|1,0,2)^, C2 : D'^> ^ 



iri,—k 



C, : |1,0,2)„ ^ -|1,0,2)^, C3 : - • (42) 

q are coefficients in an L = 2 representation, and therefore are invariant 
under Ci, (72,6*3. Ci describe continuous closed orbits in . In order to 
satisfy (^) and using the degeneracy of i?^* for k-^—k, we find that 

L = odd , k = even. (43) 

In particular, the ground state of (O) is given by L = 1, and k = 0. 



3.5 Exact Diagonalization 

The above semiclassical scheme gives a clear and intuitive picture of the 
behaviour of the system in a strong coupling limit . This limit is appropriate 
for describing, e.g., Na3 However, in Cgo the actual range of the coupling 
parameter - (? ~ 0.3 for a typical mode [p6|, E71 - suggests that the electron- 



vibron coupling is actually in the weak to intermediate regime. 

In Ref. ([0]) we have diagonalized the electron-vibron Hamiltonian (P^) 
for single Hg mode in a truncated Fock space. This approach yields accurate 
results weak to moderate values of the coupling strength. The eigenenergies 
match the semiclassical approximation ( p7|) and ( ^0|) at large g, as expected. 
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The electron-vibron pseudopotentials are defined as 



ulm = ^^+1 + K-1 - 2Er, , (44) 

where £"„ are the fully relaxed ground state energies of n electrons, which 
igbnore the effects of electron-electron interactions. For odd values of n, 
this is an effective pairing interaction often called "pair binding" in the 



literature [12 1 . In Section we found that for all odd n, the pair energies 



are negative, and given by the large g asymptotic expression 

«n=i,3,5 --g' + l-^ + 0{g-'). (45) 

The first term is the classical energy. The second term is due to reduction 
of zero point energy along the JT manifold, since only radial modes remain 
hard. This term is independent of g and positive. The last term is due to the 
quantum pseudo-rotator Hamiltonian, and the Berry phases which impose a 
finite ground state energy associated with odd L for odd numbers of electrons. 
This term, although nominally small at large g, becomes important at weaker 
coupling. If (^) is extrapolated to the weak coupling regime the last term 
would dominate the pair binding energy. The exact diagonalization indeed 
shows a significant enhancement of the pair binding energy over the classical 
value in the weak coupling regime. 

In the weak coupling limit, we can obtain analytical expressions for Un{g) 
for g « 1 by second order perturbation theory. The unperturbed Hamilto- 
nian is the non interacting part . The perturbing hamiltonian is of 
Eq. ([T^), which connects Fock states differing by one vibron occupation. All 
diagonal matrix elements vanish, and the leading order corrections to any 
degenerate multiplet are of order g"^. These are given by diagonalization of 
the matrix [E9l, 



^uLuL. = (0,n^.|i7--^^)l— iJ-|0,nL), (46) 
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in the degenerate 0-vibrons subspace. The sum imphed by the inverse op- 
erator {Ej^^ — H^)~^ extends just to the Ny=l states. The eigenvalues of 
A^^) yield the ground state energies and splittings for different electron fill- 
ings. These results, for all Hg and also Ag modes, and extended to the Ny=l 
multiplet, are discussed more extensively in Ref.|l27|. 

Here we refer only to ground state energetics. In particular, using the 
perturbative expressions, we obtain, for a single Hg, mode the small g pair 
binding energy 



flU) 



(47) 



IS 



The origin of the 5/2 factor that characterizes the perturbative result 
group theoretical [^, p^ . 

The molecular pair binding energy can be considered as an effective 
negative-U Hubbard interaction for the lattice problem, provided that the 
Fermi energy ei? is not much larger than the JT frequency scale uj. A 
mean field estimate of the transition temperature for the negative-?/ Hubbard 



model in the weak coupling regime is |12 



ep exp 



[-N{e,)\U\)-' 



(48) 



In Refs. and , the results of Migdal-Eliashberg approximation for the 
superconducting transition temperature was given. Without the Coulomb 
pseudopotentials this approach yields 



V 



wexp 



[-N{e,)\V\)-' 



(49) 



By comparing (^Vp to (^) we find a striking discrepancy between the values 
of the effective pairing interaction: 



U =3K 



(50) 
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That is to say: in the weak couphng regime, the correct molecular calculation 
yields a pairing interaction which is three times larger than the results of 
Migdal-Eliashberg theory! 



4 TDAE- Ceo: Mott-ferromagnetism 

The fullerene compound TDAEl- Cgo, where Cgo is buckminsterfullerene 
and TDAE is tetrakis(dimethylamino)ethylene C2N4(CH3)8, exhibits ferro- 
magnetism at Tc~ 16°K j^, ^. The striking aspects of this discovery are 
(i) the magnitude of T^- relatively large for a material with no transition 
metals - and (ii) its nonmetallic conductivity, suggestive of Mott-Hubbard 
localization 0. 

ESR studies ||3^ show that TDAE donates an electron to Cgo- Further- 
more, no ESR signature of TDAE+ is observed, suggesting that the TDAE 
radical spins are somehow paired. The monoclinic structure makes for a rela- 
tively short inter- Geo separation along the c-axis. We are then led to consider 
a model of CgQ chains whose conduction electrons interact via superexchange. 

Superexchange in a one-band model is always antiferromagnetic, since the 
intermediate state | t|) is a spin singlet. The molecular degeneracy of the 
tiu Geo LUMO leads to interesting possibilities not realized in a orbitally 
nondegenerate model. Indeed, Seshadri et al.^ have discussed how ferro- 
magnetism naturally arises via superexchange through intermediate states 
with a negative singlet-triplet splitting (e.g. Hund's rule) in Cgg . We intro- 
duce here what we believe to be a 'minimal model', based on the structure 
of TDAE- Cqo, which leads to insulating ferromagnetic behavior. 

The main result in this section is the full multicomponent superexchange 
Hamiltonian for the Mott insulator and the analysis of its phase diagram. 
The model is characterized by only three pseudopoentials for the intermediate 
states of Cgo in an axially symmetric crystal field. In the limit where all the 
intermediate states are degenerate there is an accidental symmetry which 
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leads to the SU(4) Heisenberg antiferromagnetic model in its fundamental 
representation. This model, coincidentally, was solved by Sutherland ||35|| 
using Bethe's Ansatz. In the case of large singlet-triplet splitting, the ground 
state is a fully polarized spin ferromagnet and an orbital antiferromagnet 
whose correlations are given by n by Bethe's wavefunction. Full details of 



the calculations below were given in Ref. [|lT 



4.1 The Hopping Model 

We consider tight binding hopping on a lattice of Cgo molecules with a filling 
of one electron per site. In general, a tetragonal or monoclinic crystalline 
symmetry will resolve the triply degenerate tiu orbital into three distinct 
levels. Better details could be obtained ab-initio once the precise structure 
and orientations of the Cgo molecules are experimentally ascertained. In our 
model, we shall retain only what we believe may be the essential microscopic 
physics undelying the ferromagnetism in TDAE- Cqq. (a) The hopping is 
quasi-one dimensional along the c-axis. (b) We assume that the crystal field 
resolves the tiu orbital triplet into a lower doublet (/ = ±) and a higher singlet 
/ = at higher energy, as if the crystal fields are cylindrically symmetric 
about an axis which pierces the center a pentagonal face of Cgo- (c) Hopping 
along the chains is assumed to preserve the orbital magnetization I. 
Thus we investigate the Hamiltonian Ti = Ti^op + Ti-^^p + '^ion, where 

hop — X! ('^^^'^ 

^hop = -iE4(^±)(c^''^Wqv(i + '5±) + H.c. 

Hion = E^A |A(i)) (A(i)l . (51) 

Here c^'°'(i) creates, at site i, an electron of spin polarization a =t, | and 
"isospin" / = -|-, — . c and S± denote nearest neighbor lattice vectors in the c 



^hop = -t\\j:{c^'^ii)ci,ii+c)+ii.c. 
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direction and the a-b plane respectively, tu/ are the hopping matrix elements 
between orbitals / and I' on neighboring chains. 

Tiion is the interaction Hamiltonian which discourages multiple electron 
occupancy on any Ceo niolecule. It is parametrized by pseudopotentials u\ 
which correspond to the following Cgg multiplets, 

u^: _i.(ct+Tct-i_ct+ict-T)|0) 

ui : c^+^c^-^ I ) , ^(c^+^ct-i + c^+^ct-T) I ) , c^+^c^-^ \ ) 

U2 : c"^+Tc"^+^ I ) , c"^~Tc"^~^ I ) (52) 

The relations between ma and the isotropic pseudopotentials ul of angular 
momenta L are: Ui = Ui, U2 = U2, but uq = |mo + |m2 due to projecting out 
the "/ = 0" orbital state. Thus, while in an isotropic environment there might 
be pair binding {uq < 0) due to electron-electron [|l^] and electron-vibron [|10| 



interactions, it does not preclude a repulsive > in the monoclinic crystal 
field environment. This may help to explain why TDAE- Ceo is not a CDW, 
nor a superconductor as is AaCeo- 



4.2 Multicomponent Superexchange Hamiltonian 

Experiments have shown that TDAE- Ceo is insulating at low tempera- 
tures, consistent with the existence of a gap to charge fiuctuations (i.e. all 
Ma > 0) 0. The low- lying excitations can be described by a superexchange 
Hamiltonian, formally obtained as a second order expansion in small t"/M. 
Since charge excitations are gapped, a renormalized version of the superex- 
change Hamiltonian is expected to describe the low energy excitations also 
for tyu>l. 

The zeroth order states of the superexchange Hamiltonian are four singly 
occupied states enumerated by /, a. The operators which act on these states 
can be represented by spin operators, = \Y.i,a,u' f^'^^" {i) 'T'aa' ^la'i^) ^'^d 
"isospin" operators. If = \ J2i,i',ct c^''^(i) tI'i, Ci,^{i) where r are the Pauli ma- 



19 



trices. Taking into account the constraint J2a c^'^c^ = 1, the 15 independent 
elements of the SU(4) generators 5*^ = c^"c^ can be expressed in terms of 
the 15 operators {Sf", T , Sf'F}. 

For simphcity we consider the purely one-dimensional limit where t^, = 0. 
There are three superexchange constants defined as 

J^ = ^^ M = 0,1,2, (53) 

Um 

A straigthforward, though cumbersome, leads to an effective Hamiltonian 
given by 



where 



H = ■ Sn+l + B In ■ In+l + C rj^^^+ (54) 

n 

+D Sn ■ Sn+1 In ' In+1 + E Sn ■ Sn+1 I^In+l + -^y 
A = —\Jl + J2 + \ J0) B = I Jl — i Jq; 

C = Jo - J2, D = 2Ji + 2Jo, 

E = AJ2- 4Jo, F = -| Jl - i J2 - I Jo . (55) 

This model possesses a global SU(2)xU(l) symmetry, i.e. H commutes with 
J2n and with J2n ^n- Enlarged symmetries occur when Jo = J2, where 
the symmetry group is SU(2) xSU(2), and when Jo = Ji = J2, where the 
symmetry group is SU(4). 

SU(4) Point - At the point ui = U2 = uq = u (^) acquires full SU(4) 
symmetry. For each c-chain the Hamiltonian is 

H5C/(4) = JEE^|W^f(^+l)' (56) 

where J = 2{t^^/u. 

The SU(P) Heisenberg antiferromagnet in the fundamental representa- 
tion has been solved by Sutherland for general P using Bethe's Ansatz 
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This model exhibits P — 1 gapless elementary excitation branches. We pre- 



sume, based on what happens in the SU(2) model |^6[, that for a chain of 
sites where N is an integer multiple of P, the ground state is an SU(P) 
singlet and the low-lying excitations transform according either to the singlet 
or the adjoint representation. This is essentially what happens in the fermion 
mean field (large P) theory of the SU(P) antiferromagnet |3^, The 
mean field has four degenerate quarter-filled (fcp = ^vr) bands for P = 4. 
Although there it has no true long-ranged order, the spin and isospin sus- 
ceptibilities diverges at the nesting wavevector 2kp = ^vr, which describes a 
commensurate spin density wave of period four. The period four arises be- 
cause the spin chain is in its fundamental representation, and by '4-ality' one 
needs four sites to make a singlet [^. The mean field theory also predicts 
a constant uniform (Pauli) susceptibility, and a linear specific heat as in a 
Fermi hquid[^ . 



Ferro-Antiferromagnetic points - Along the surface U2 = Uq, our Hamilto- 
nian possesses an SU(2) xSU(2) symmetry. There are then two special limits 
in which we can determine the exact ground state, (i) The "F x A model" 
at Mo — ^ oo, with = 2(tll)^/-ui, 

^FxA = -^^Y^i^n ■ Sn+1 + " In " ^n+l) (57) 

Ml „ 

where the interactions are ferromagnetic in the spin channel and antiferro- 
magnetic in the isospin channel, and (ii) the "A x F model" for Mi — oo, 
with J|| = 2(tl')^/'Uo, and the roles of I and S interchanged. 

It is possible to prove that the ground state of TYfxA is the fully polarized 
ferromagnet | F ) ^ for the spin variables, and Bethe's ground of the spin-half 
antiferromagnet for the isospin variables i.e. 

^FxA = \ F)s® iBethe)^ . (58) 

A corresponding result holds for TYaxFj with spin and isospin variables ex- 
changed. Due to the SU{2) (g) SU{2) symmetry the total spin S'tot, total 
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isospin Jtot, and their polarizations along the z axis {Ms and M/, respec- 
tively) are good quantum numbers. Following Lieb and Mattis' proof of 
the Marshall theorem for the Heisenberg model PT|, we perform a vr rotation 



about the z axis of the isospin operators on odd-numbered sites. The Hamil- 
tonian transforms into a non-positive ('negative semidefinite') operator in 
the product Ising basis 

n 

x{Sn-Sn+l + l)=n'^x^ (59) 

The accessibilty of all states within a given magnetization sector by repeated 



application of the Hamiltonian, implies (see Ref. |^) that the ground state 
of HpxA in the sector (M5, M/) = (0, 0) can be chosen to be positive definite 
in the sublattice-rotated Ising basis, i.e. it obeys Marshall's sign rule. Since 
the same Marshall signs hold for the state on the right hand side of Eq. 
|58| , which has Stot = and Jtot = 0, the two sides of Eq. ^ have finite 
overlap hence the same ^tot and Jtot- We are free to choose Ms = as a 
representative of the ground state manifold. Note that | \E'o ) is indeed an 
eigenstate of the spin triplet projection operator (Sn-Sn+i + j) with eigenvalue 
one. It follows from Eq. ^ that the isospin part of the wavefunction is the 
ground state of the spin-half antiferromagnetic Heisenberg chain, given by 
Bethe's Ansatz. 

Exact excitations of Tipx a within the isospin sector (retaining full spin po- 
larization) with dispersion |vr J| sin k\ can be constructed as in Refs. ||3^, ^ 



The gapless ferromagnetic magnons, which exist due to Goldstone's theo- 
rem, can be approximated within the Single Mode Approximation (SMA): 
I A; ) = 5^ I \I'o ) . The trial state dispersion is 

uj{k) < 2 ln(2) J(l - cos k) (60) 

from which we see that the ferromagnon bandwidth is decreased due to the 
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antiferromagnetic nearest-neighbor isospin correlations, i.e. (| — I„ ■ In+i) = 
ln(2). 

4.3 Classical Phase Diagram 

The ground state depends on the dimensionless ratios Uq/ui and U2/U1. The 
classical approximation (justified at S,I » 1) is given by minimizing the 
bond energies of Eq. (|5^ ) as function of vectors Si and Jj of magnitude |. 
The results are plotted in Fig. |^). 

It is interesting to note that the SU(4) symmetry point is at the border 
of 4 distinct ordered phases of different symmetries, where the energy is 
degenerate along the lines (I„, ■ In+i) = — | and (S„ ■ S„+i) = |. The large 
degeneracy of the classical SU(4) model is reduced by quantum fluctuations. 

The classical regime of Heisenberg spin-ferromagnetism and isospin-anti- 
ferromagnetism extends throughout uq,U2 > Ui, although quantum fluctua- 
tions break the SU(2) isospin symmetry away from the isotropic line Uq = U2 
(marked dashed line in Fig. H). 

3D Ordering in the F x A Model - As shown by Scalapino et al. |[43|| , 



one 



can treat the interchain interactions by mean fleld theory and thereby derive 
an expression for the full susceptibility Xab^lLiIzi^) in terms of x]^{qz,oj), 
the susceptibility for the one-dimensional chains. The general result is 

X(g±, g„ oj) = [l- J±{q^) X^°(g„ u)] x'^(g„ u) , (61) 

where J±{q±) = J2d^ J±i^±) e~'^'^^'^^ is the spatial Fourier transform of the 
interchain coupling matrix. (Note that the quantities X, J_l, and in 
Eq. ^ are matrices.) This approximation also may be employed at flnite 
temperature. 

Consider now the Fx A model discussed above. At flnite temperature T, 
long-ranged ferromagnetic order is destroyed and the global SU(2)xSU(2) 
symmetry is restored. The uniform susceptibility of the ferromagnetic chain 
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is given by 

^f(0'0;^) = ^ + -- - ' (62) 



as was first derived by Takahashi in Ref. (see also Refs. ||38|, 0). 

For the antiferromagnetic susceptibility, we appeal to the bosonization 
results of Schulz and of Eggert and Affleck who have computed the 
dynamic susceptibility of the S = ^ antiferromagnetic Heisenberg chain. 
Performing a Fourier transform of their result and taking the low frequency 
and wavevector limit near the antiferromagnetic point we obtain the staggered 
isospin susceptibility 

X. « ^ (63) 

where Oq — 4.44. 

For mixed interchain coupling operators e.g. O = S^P, we may use the 
assumed independence of ferromagnetic and antiferromagnetic magnons to 
obtain at low temperatures Xpj^{7r,0;T) ~ (JyT)"^/^, which diverges even 
more slowly than Xa in the T —* limit. 

The interchain interaction is given by J_l = J||(()f:]'-)^ + (t^)^)/4(tll)^, where 
t^,t2 are the transverse hopping integrals (see Ref. [0]). Thus, as the tem- 
perature is lowered, a transition from paramagnetic to ferromagnetic state 
should set in when J±X'p = 1. This yields Tq ~ JJhJ±/24. The relation 



Tc oc J\\J± was also found by Scalapino et al. (ref. fi^) in their studies of 
anisotropic Heisenberg magnets. It is conceivable that at still lower temper- 
atures a Neel ordering of the isospin variables occurs at a Neel temperature 

4.4 Experimental notes 

(a) The lower isospin transition, to our knowledge, has not been resolved 
experimentally. Perhaps it is not very well separated from the ferromag- 
netic transition which would help explain the mysterious excessive entropy 
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of transition found by Ref.p^ 



(b) Alternatively, the isospin ordering might be preempted by a isospin- 
Peierls ordering (orbital dimerization) aided by the electron-phonon cou- 
pling. In that gnature for the isospin-Peierls effect should be present 
in X-ray scattering or in the phonon spectrum. 

(c) The role of a possible orientational disordered ground state has 
not been considered here although it might help explain the observed weak 
ferromagnetism 0. In addition, Bloch's T^/^ temperature dependence of 
the ordered moment found in Ref. [Q which holds upto T Tc is hard to 
reconcile with quasi-one dimensionality where J± << Tc- 

In summary, this section has described a model of quasi-one dimensional 
interacting electrons with doubly degenerate orbitals motivated by the struc- 
ture of TDAE- Ceo- At occupancy of one electron per site, we obtain a Mott- 
insulator with muticomponent superexchange between spins and isospins at 
neighboring sites. At special values of the interactions we identify exactly 
solvable points, including the SU(4) antiferromagnet, and spin-ferromagnet, 
isospin-antiferromagnet limit. The classical ground state diagram also con- 
tains a large region of spin ferromagnetism and orbital antiferromagnetism 
which we believe is relevant for TDAE- Cgo- A mean field analysis of the 
interchain coupling in this regime predicts two transition temperatures: fer- 



romagnetic spin ordering at Tq oc ^J\\J±, and orbital (isospin) antiferromag- 
netic ordering at c<. J±. This lower transition, to our knowledge, has not 
yet been resolved experimentally. 



5 Concluding remark 

In this review we have seen that it is possible to understand some of the 
unusual electronic properties of buckminsterfuUerene compounds by relating 
them to the high molecular symmetry, and local degeneracies of the partially 
filled tiu orbitals. We have investigated in detail some simplified models with 
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local degeneracies, which exhibit enhanced superconductivity, and ferromag- 
netism, as well as other possibile phases not yet observed experimentally. 
This suggests that continual experimenting with the family of fuUerene com- 
pounds would most probably produce further surprises. 
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Table 1: Semiclassical ground state distortions and energies for a single Hg 
coupled mode of frequency u. n is the electron number, S is the total spin, 
Zn,fn are the JT distortions, Ui is the occupation of orbital i, En is the 
ground state energy and, and m„ is the pair energy (Eq.(^^). Energies are 
calculated for strong coupling to order g'"^. 
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Figure 1: A polar representation of the Jahn- Teller distortions w^"^ {6 , (p) , 
Eq. (P^) ■ The distortion is measured relative to a sphere, (a) The unimodal 
distortion for the ground states of n = 1,2,4,5 electrons (b) The bimodal 
distortion for n = 3 electrons. 



Figure 2: TDAE- Cgo: Classical ground state phase diagram 
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